Libraries:
library(tidyverse)
library(readr)
library(rvest)
library(stats)
library(readxl)
library(dplyr)
library(stringr)
library(ggplot2)
library(ggthemes)
library(stringr)
library(data.table)
library(geojsonio)
library(leaflet)
library(rgdal)
library(haven)
library(stargazer)Data:
Model 3:
First making a graph of annual GDP highs and lows – maybe a temporary proxy for recessions?
Annual enrollment graph:
Join enrollment data and gdp data to create linear model test:
Graph it?
College Proximity Question 5/3: (Reading in Ivy’s data)
cz <- read_dta('cz.dta')
Error in read_dta("cz.dta") : could not find function "read_dta"Read in/create mobility data: (Trends in Mobility: Commuting Zone Intergenerational Mobility Estimates by Birth Cohort) https://opportunityinsights.org/data/?geographic_level=101&topic=0&paper_id=0#resource-listing
Read in geojson file:
Commuting zones on the map (cz.geojson@data) are in 1990s format. They need to be converted so our post-2000 data can be connected to the shapefiles: (https://www.ers.usda.gov/data-products/commuting-zones-and-labor-market-areas/)
Making the map actually display data:
cz.geo %>%
leaflet() %>%
addTiles() %>%
addPolygons(fillColor = ~colors1(cz.geo@data$ncollege),
weight = 1,
color = "white",
opacity = 0.5,
fillOpacity = .7) %>%
setView(-96, 37.8, 3) %>%
addLegend(pal = colors1,
values = cz.geo@data$ncollege,
title = "Number of Colleges in Commuting Zone")
Some values were outside the color scale and will be treated as NA